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Abstract 

In this paper we report the first measurement of e^e~ £)(*)+£)(*)- processes. The cross- 
sections for e+e~ D*^D*~ and e+e~ D^D*~ at -y/s = 10.58GeV/c^ have been measured to 
be 0.65 ± 0.04 it 0.07pb and 0.71 ± 0.05 ± 0.09 pb, respectively. We set an upper hmit on the cross- 
section of e"'"e~ D'^D^ of 0.04 pb at tlic 90% confidence level. In addition we have measured the 
fraction of the Dl^D*^^ final state in the e+e" D*+D*- reaction to be (97 ± 5)%. The analysis 
is performed using 88.9 fb~^ of data collected by the Belle detector at the e+e~ asymmetric collider 
KEKB. 

PACS numbers: 13.65.+i, 12.39.Hg, 13.87.Fh 
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The processes e+e" D^*^D^*\ with no extra fragmentation particles in the final state, have 
not previously been observed at energies a/s 3> 2M£,. The cross-sections for these processes 
can be derived once the charmed meson form factors are determined for the appropriate 
value of momentum transfer, = s. In the HQET approach based on the heavy-quark spin 
symmetry, the heavy meson form factors are represented in terms of a universal form factor, 
called the Isgur-Wise function. However, for asymptotically large 3> M^/A the leading- 
twist contribution, which violates the heavy-quark spin symmetry, becomes dominant 
For an intermediate range of momentum transfer, sub-leading twist corrections are also 
important. A calculation that takes these effects into account (H predicts cross-sections 
of about 5pb"^ for e+e" DD* and e+e~ — > D^Dl at ^ ~ 10.6 GeV (the subscripts 
indicate transverse [T] and longitudinal [L] polarization of the D*); the cross-section for 
e~^e~ — > DD is expected to be suppressed by a factor of ~ 10~^. 

In this paper, we present the first observation of the high energy e~^e^ D*^D*^ and 
e^e^ —>■ D^D*^ processes, and a measurement of their cross-sections and polarizations. We 
also set an upper limit on the cross-section for e~^e~ —>■ D^D~ . 

The present study is limited to final states that contain charged D*^*-* mesons only: this 
simplifies the analysis, since it is possible to select events with no neutral particles in the 
final state. Since the contribution of electromagnetic current coupled to light quarks is 
negligible compared to that for heavy quarks, the neutral and charged charm cross-sections 
are expected to be the same p*]. 

The analysis is based on 88.9 fb^^ of data at the T(4S') resonance and nearby continuum, 
collected with the Belle detector [3] at the KEKB asymmetric energy storage rings i. We 
select well-reconstructed tracks consistent with originating from the interaction region as 
charged pion candidates. Those passing particle identification cuts based on dE/dx, aerogel 
Cerenkov, and time-of-fiight information are selected as charged kaon candidates. We 
then reconstruct D° and D~^ mesons in the decay modes D^ — > K~7t~^, D^ K'h^tt^tt' 
and D^ K^-k^-k^. The selected combinations are constrained to a common vertex, 
and quality cuts are imposed on the vertex fit to reduce the combinatorial background. 
A 15MeV/c^ interval around the nominal D masses is used to select D^ K~7r~^ and 
D~^ K^TT^Tc'^ candidates; for the —>■ K~tt^'k^tt~ decay mode the signal window is 
chosen to be lOMeV/c^ around the nominal D^ mass. The selected D candidates are then 
subjected to a mass and vertex constrained fit to improve their momentum resolution. The 
D*^ mesons are reconstructed in the D^ir^ decay mode. The invariant mass of the D^t:^ 
combination is required to be within a 2MeV/c^ (~ 3 a) mass interval around the nominal 
D*^ mass. 

The processes e"*"e^ [){*)+[)(*)- can be identified by energy-momentum balance in 
fully reconstructed events that contain only a pair of charm mesons. However, the small 
charm meson reconstruction efficiency via the studied channels results in a very small total 
efficiency in this case. Taking into account two body kinematics, it is sufficient to reconstruct 
only one of the two charmed mesons in the event to identify the processes of interest. For 
simplicity, we refer to the fully reconstructed D meson as the D^*^~^, and the other as the 
D*~; the charge-conjugate modes are included in the analysis. We choose the mass of the 
system recoiling against the reconstructed (MrecoiK-D'-*-'^)) as a discriminating variable: 

where -y/i is the total center of mass (CMS) energy, and Ej^(,)+ and 'Po(*)+ are the CMS 
energy and momentum of the reconstructed D^*^~^. For the signal a peak in the Mrecou 
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distribution around the nominal mass of the recoihng D~ or D*^ is expected. This method 
provides a significantly increased efficiency, but also a higher background, in comparison to 
full event reconstruction. For the e^e~ — > D~^D*~ and e~^e~ — > D*~^D*~ processes we find a 
better compromise between higher statistics and smaller background: the first is fully 

reconstructed, while the recoiling D*~ is required to decay into D^t^^iow- The reconstructed 
'^dow provides extra information that allows us to reduce the background to a negligible level 
as explained below. 

We calculate the difference between the masses of the systems recoihng against the 
£)(*)+7r^o^ combination, and against the alone. 

The variable AMrecoii peaks around the nominal D*^ — mass difference with a resolution 
of (TAMrocoii ~ IMeV/c^ as found by Monte Carlo simulation. For e+e~ D*^D*~ and 
e+e" — > D~^D*~we combine D*~^ and candidates together with tt^q^^ and require AMrecoii 
to be within a ±2 MeV/c^ interval around the nominal Md*+ — M^o mass difference. 

The Mrecoii(-D*''') and Mrecoii(-D''") distributions are shown in Figs, la and lb, respectively. 
Clear signals arc seen around the nominal D*~ mass in both cases. The higher recoil mass 
tails in the signal distribution arc due to initial state radiation (ISR). The hatched histograms 
show the Mrecoii distributions for events in the AMj-ecoii sidebands. 

The backgrounds in the region of -Mrecoii < 2.1 GeV/c^, are negligible for both processes. 
To provide a numerical estimate, we divide the background sources into three categories: 

I fake reconstructed D*'^ or 

n-K, where the vr^^^^ is not produced from D* decays, and, thus, 
produces no peak in the AMiccoii distribution; 

III e+e" D'^*^+D*-mr, where n > 1. 

First we consider the backgrounds for the process e+e" D*^D*~ . To estimate back- 
ground (I) we count the entries in the Mrecoii (-D^tt"*") < 2.1 GeV/c^ interval for D^vr"*" combi- 
nations taken from the D*^ mass sideband (2.016 < M£)Q^+ < 2.02GcV/c^). Three events 
are found in the data, while the Monte Carlo predicts a contribution of 2.5 events the from 
non-Gaussian tails in the Mr)*+ resolution function. The signal Monte Carlo is normahzed 
to the number of entries in M^ccoii{D*~^) < 2.1 GeV/c^ region in the data. Background (II) 
is estimated using AMrecoii sidebands (0.150 < AMrecoii < 0.154 MeV/c^). In the region 
-Mrecoii (-D*''') < 2.1GeV/c^ 8 events are found in the data; according to Monte Carlo, 4 
events are expected from initial state radiation. Thus, backgrounds (I) and (II) are esti- 
mated to be smaller than 3 and 9 events at 90% CL, respectively. As a cross-check we also 
study the wrong-sign D^*'''^ tt^iow combinations: 2 events with wrong-sign tt^^^'s are found 
in the data in the interval Mrecoii(-D*+) < 2.1 GeV/c^. 

The remaining background (III) can result in peaks in both the M(£)*+) and AMrecoii dis- 
tributions, but has a threshold in Mrecoii distribution at Mi,*+ + M^o = 2.15GeV/c^, which 
is ~ 1(7 away from the chosen MrccoiK-D*"*") interval. To estimate the residual background 
(III) contribution we perform a fft to the Mrecoii (-D*"*") distribution. The signal function is 
the result of convolving the generated Mrecoii(-D*"*') distribution with the detector resolution: 
the signal function is the sum of a core Gaussian and an asymmetric function represent- 
ing the Mrecoii ("D*"^) shape when the studied process is accompanied by radiative photon(s) 
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FIG. 1: a-b) The distributions of the mass of the system recoihng against a) D*~^ , and b) . 
Points with error bars show the signal AMrecoil region; hatched histograms correspond to AMrecoii 
sidebands. The sohd hnes represent the fits described in the text; the dashed Unes show the 
contribution due to events with ISR photons of significant energy. The dotted hnes show the 
expected background contribution, c) The distributions of Mrecoii(^^) without any requirement 

on AMrecoil- 

with significant energy. The M^ecoi\{D*'^) resolution due to detector smearing and the signal 
function offset are left as free parameters in the fit, to test the agreement with the Monte 
Carlo predictions. The background (III) distribution is parameterized by a threshold func- 
tion, a ■ (Mrecoii(-D*+) — M{D*^)pDG — M (7r°)p£)G)^5 convolved with the detector resolution, 
where a and (3 are free parameters. The fit results are shown in Fig. ^p, as the solid line. 
The dotted line represents the expected background (III) distribution. The signal yield is 
found to be 815 ib 28 events. The Mrecoii resolution a = 56.1 ± 2.2MeV/c^ is found to be 
in excellent agreement with the Monte Carlo expectation (56.4MeV/c^), and the shift of 
the signal peak position in the data with respect to the Monte Carlo position is found to 
be consistent with zero (0.6 ± 2.5MeV/c^). The contribution from background (III) in the 
-^recoil < 2.1 GcV/c^ interval is estimated from this fit to be less than 2 events. 

The backgrounds for the e~^e~ D^D*^ process are studied in a similar way. Five 
events are found in the data in the mass sideband region (20 < \MK-Tr+n+ — M^+l < 
35MeV/c^), whereas 1.5 events are expected from the signal Monte Carlo. There are 6 
data events in the AMrecoii sidebands with Mrecoii(-D^) < 2.1 GeV/c^, while 6 signal events 
are expected from the Monte Carlo simulation. Thus backgrounds (I) and (II) for the 
e'^e^ D^D*^ process are estimated to be smaller than 7 and 4 events at 90% CL, 
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respectively. Three wrong sign D^nfi^^ combinations are found in the data in the interval 
Mj-ecoi\{D~^) < 2.1GeV/c^. A similar fit is then performed to the Mj-ecoi\{D~^) distribution. 
The signal yield is found to be 423 ± 20 events. Again, excellent agreement between the 
Monte Carlo and signal shape parameters is found: a = 58.1 ± 3.6MeV/c^ (60MeV/c^ 
is expected from the Monte Carlo); M{D*-)data - M{D*-)mc = -2.1 ± S.eMeV/c^. We 
conclude that background (III) is smaller than 2 events. 

From the above study we estimate the total background in the Mrecoii < 2.1GeV/c^ 
interval to be smaller than 14 and 16 events for the e^e~ —>■ D*~^D*~ and e^e^ D^D*~ 
processes, respectively, which is of order of 1% of the signal. We assume all events in the 
interval Mrecoii < 2.1 GeV/c^ are signal and include the possible background contribution in 
the systematic error. 

Since the reconstruction efficiency depends on the production and D*^ helicity angles, 
we perform an angular analysis before computing cross-sections. The helicity angle of the 
non-reconstructed D*~ is calculated assuming two-body kinematics. A scatter plot of the 
hehcity angles for the two D*-mesons from e+e" D*^D*~ {cos (f){D*^^) vs cos0(D*Q„_^g^)) 
for the recoil mass region Mrecoii (-D*''') < 2.1GeV/c^ is shown in Fig. I^b. The distribution 
is fitted by a sum of three functions corresponding to the D^D^, D^D^ and Dj^D^ final 
states, obtained from Monte Carlo simulation. The fit finds the fractions of D^D^, D^DI 
and DlDl final states to be (1.5 ± 3.6)%, (97.2 ± 4.8)% and (1.3 ± 4.7)%, respectively. 
Figure shows the D*~ meson helicity distribution for e~^e~ D^D*~. The fraction of 
the D'^D^ [6J final state is found from the fit to be equal to (95.8 ± 5.6)%. 
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FIG. 2: a) The scatter plot cos{(j)D*^^) vs cos(</>£)*^^ ^^^) {e^e — > D*^D* ). b) D*^ meson helicity 
angle distribution for (e'^e" D^D*^) signal candidates. 

The raw production angle distributions for D*^ from e+e^ — » D*^D*^ and from 
e^e~ — > D~^D*~ processes are obtained from the region of recoil masses Mrecoii < 2.1 GeV/c^. 
In this region, the efficiency for signal events with initial state radiation photon(s) with sig- 
nificant energy is low. We calculate the fraction of rejected ISR events according to a Monte 
Carlo simulation based on Ref. 0], in bins of the production angle. The reconstruction 
efficiency is estimated from the Monte Carlo simulation, in the same production angle bins, 
taking the measured D*"^ polarization into account. Figures 3a and 3b show the production 
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angle distributions for D*~^ from e~^e~ D*~^D*~ and from e~^e~ D~^D*~, respec- 
tively, after correcting for the reconstruction efficiencies and intermediate branching ratios. 
These distributions are fitted with the function iV • (1 + acos^^). To calculate the total 
cross-sections, the signal yields are corrected to take into account the fraction of events with 
initial state radiation that lie outside of the interval Mj-ccon < 2.1 GeV/c^. The efficiency cor- 
rected signal yields are found to be 58000 ± 3400 and 64000 ±4800 for the e+e" D*+D*- 
and e^e~ D~^D*~ processes respectively. The parameters a are equal to 0.8 ± 0.3 and 
2.3lo;7 for the two processes. 
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FIG. 3: The production angle distribution for a) e+e~ D*~^D*~ and b) e+e~ D^D*~ 

processes after correcting for reconstruction efficiencies and intermediate branching fractions. 

We find cross-sections of 0.65 ± 0.04 ± 0.07 pb and 0.71 ± 0.05 ± 0.09 pb for e+e" 
D*+D*- and e~^e~ D^D*~ , respectively, where the first error is statistical and the second 
systematic. The sources of systematic error are summarized in Table EJ The dominant 
contributions are from the uncertainties in tracking efficiency and corrections for the initial 
state radiation. 



Source 


e+e- ^ D*+D*- 


e+e- D+D*- 


Tracking efficiency 


9% 


8% 


Identification 


2% 


2% 


Backgrounds 


n% 


1J% 


ISR correction 


5% 


5% 




4% 


8% 


Total 


11% 


13% 



TABLE I: Contributions to the systematic error for the e^e D*^D* and e^e D^D* 
processes. 

We search for the process e^e" D^D^ by studying the recoiling against the recon- 
struted {Mrecou)- In the e"'"e~ D*^D*^ and e^e~ D^D*~ analyses, backgrounds 
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are strongly suppressed by the tight AMrecoii cut, which is not apphcable dor e+e~ D~^D~. 
Without this cut the combinatorial background is significant (~ 20%). We use mass 
sidebands (20 < IMk-kit — Md\ < 35MeV/c^) to extract the Mrccoii distribution for the 
combinatorial background. Fig. Ic shows the distribution of MrecoiK-D^) after sideband sub- 
traction. To extract the e~^e^ D^D~ and e+e^ — D^D*^ yields we fit this distribution 
with the sum of two signal functions corresponding to D~ and D*" peaks and a background 
function. The latter is a threshold function, a ■ {x — M{D~)p£,g — M{7i^)pdg)'^ , convolved 
with the detector resolution, where a and f3 are free parameters. For the fit we use only 
the region Mrecou < 2.25GeV/c^, because of the contribution of e'^e" [)(*)[)** at higher 

-^^recoil • 

The fit finds — 13±24 events in the peak and 935±42 in the D*'^ peak. The fit function 
is shown in the Fig. Ic as the solid line; the dashed line shows the contribution of events 
with ISR photons of signficant energy (larger in this case due to absence of the AMrecoii 
cut); and the dotted line represents the case where the contribution of e+e^ D^D~ is 
set at the value corresponding to the 90% confidence level upper limit. The reconstruction 
efficiencies for the e~^e~ — > D^D~ and e+e~ — * D~^D*~ are found from Monte Carlo to be 
3.1% and 1.7%, respectively. In the e~^e^ D^D^ Monte Carlo the production angle is 
assumed to have a sin^ 9 distribution, as required by conservation of angular momentum; for 
e~^e~ D^D*~ the production angle distribution is fixed according to the analysis presented 
above. The e'^e" — >• D^D*~ cross-section is calculated to be 0.61 ± 0.05 pb which agrees 
with the result using the AMrecoii method. However, in this case the systematic uncertainty 
in the signal yield is large due to the non-negligible e^e~ — D^Dtt background under the 
peak, which can only be determined from the higher Mrecoii region. For the e"'"e" D^D~ 
cross-section we set an upper limit of 0.04 pb at the 90% confidence level. 

In summary, we report the first measurement of the cross-sections for the e~^e^ — > D*^D*~ 
and e+e" ^ D+D*- processes at ^/i = lO.GGeV to be 0.65±0.04±0.07pb and 0.71±0.05± 
0.09 pb, respectively, and set an upper limit on the e+e^ D^D^ cross-section of 0.04 pb 
at 90% confidence level. The measured cross-sections are an order of magnitude lower than 
those predicted in the Ref. jH, but their relative sizes are as predicted: the cross-sections 
for e^e~ D*~^D*~ and e^e~ D^D*~ are found to be close each other, while the cross- 
section for e+e^ D^D^ is much smaller. The helicity decomposition for e~^e^ D*^D*~ 
is found to be saturated by the D'^D*]^ final state (the fraction is equal to (97.2 ± 4.8)%) 
and for e+e~ — >■ D^D*~ - by the final state (95.8 ± 5.6%), in good agreement with the 
predictions of Ref. jH . 
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